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Theoretical Interpretation of the Measurements of the Secondary 
Eclipses of TrES-1 and HD209458b 

A. Burrows^, I. Hubeny^, & D. Sudarsky^ 

ABSTRACT 

We calculate the planet-star flux-density ratios as a function of wavelength 
from 0.5 /imto 25 /imfor the transiting extrasolar giant planets TrES-1 and 
HD209458b and compare them with the recent Spitzer/IRAC-MIPS secondary 
eclipse data in the 4.5, 8.0, and 24 /rm bands. With only three data points and 
generic calibration issues, detailed conclusions are difficult, but inferences regard¬ 
ing atmospheric composition, temperature, and global circulation can be made. 

Our models reproduce the observations reasonably well, but not perfectly, and 
we speculate on the theoretical consequences of variations around our baseline 
models. One preliminary conclusion is that we may be seeing in the data in¬ 
dications that the day side of a close-in extrasolar giant planet is brighter in 
the mid-infrared than its night side, unlike Jupiter and Saturn. This correspon¬ 
dence will be further tested when the data anticipated in other Spitzer bands are 
acquired, and we make predictions for what those data may show. 

Subject headings: stars: individual (TrES-1, HD209458)—(stars:) planetary 
systems—planets and satellites: general 


1. Introduction 

High-precision radial-velocity measurements of many nearby stellar primaries have re¬ 
vealed the presence of extrasolar giant planets (EGPs) in orbit around them (Mayor & Queloz 
1995; Marcy & Butler 1996; Marcy & Butler 1998; Marcy, Cochran, & Major 2000, and ref¬ 
erences therein) From these data, the projected mass (MpSin(z), where Mp is the planet 
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mass and i is the orbital inclination) and orbital parameters of more than 150 planetary 
companions have been determined. These data have revolutionized the study of planetary 
systems, not the least because their orbits and masses vary widely and are in general quite 
unlike those of our Jovian and ice giants. 

The detection of a few transiting EGPs (currently seven), whose orbits are nearly edge- 
on, has provided radii, masses, and inclinations (Henry et ah 2000; Charbonneau et al. 2000; 
Brown et al. 2001; Sasselov 2003; Konacki et al. 2003ab; Konacki et al. 2004; Bouchy et al. 2004; 
Pont et al. 2004; Torres et al.2005). For HD209458b, there are in addition indications of 
its atmospheric composition from the wavelength-dependence of the photometric dip of the 
stellar light during transit (Charbonneau et al. 2002; Vidal-Madjar et al. 2003; §2). Figure 1 
depicts the mass-radius plot for these transiting EGPs, with Jupiter and Saturn themselves 
superposed. 

However, it is only by direct detection of the planets and their photometric and spec¬ 
troscopic characterization that they can be studied in depth to reveal their atmospheric 
and physical properties (Sudarsky, Burrows, & Hubeny 2003; Allard et al. 2003; Burrows, 
Sudarsky, & Hubeny 2004; Burrows 2005). While orbital distances (a), periods (P), and 
eccentricities (e) are reasonably well determined, the investigation of an EGP in physical 
detail requires at a minimum the actual mass (Mp), radius (Rp), and composition. Space 
provides the needed access, but it had generally been thought that high-contrast imaging 
is necessary to separate out the planet from the bright star. However, Charbonneau et al. 
(2005) and Deming et al. (2005) have recently shown that variations in the summed light 
of the planet and star for the close-in transiting EGPs TrES-1 and HD209458b can be de¬ 
tected with Spitzer during secondary eclipse. Secondary eclipse is ~180° out of phase with 
the transit and is when the planet is occulted by the star, thereby shutting off the planet’s 
contribution to the summed light. The approximate magnitude of this diminution varies 
signihcantly with wavelength, with the mid- to far-infrared being the most favorable bands. 
At 10-30 microns and in the Spitzer/IRAC bands from 3.6 /rmto 8.0 jumthis contrast was 
predicted to be ~ 10“^ (Burrows, Sudarsky, & Hubeny 2003; Burrows 2005). This is near 
what has now been detected. In this paper, we customize to the TrES-1 and HD209458 
systems the calculation of the planet-star flux-density ratio as a function of wavelength and 
compare the resulting theory with the Spitzer secondary eclipse data to draw conclusions 
about close-in EGP atmospheres. In §2, we summarize the new Spitzer data and in §3 we 
briefly describe our numerical techniques. Then, in §4 we present our results from 1 /rmto 
25 /imand wrap up in §5 with a discussion of our conclusions and the outstanding issues 
concerning irradiated EGPs that such eclipse data might address. 
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2. Summary of TrES-1 and HD209458b Data 

The parameters of a transiting planetary system determine the inpnts to a theoretical 
calcnlation of its spectrnm and its planet-star flnx-density ratio as a fnnction of wavelength. 
Incident radiation at the planet’s snrface is a fnnction of the stellar flux density and the 
orbital distance, and the planet’s mass and radius determine the gravity of its atmosphere. 
For both TrES-1 and HD209458b, the interior flux (Burrows et al. 2000) is dwarfed by the 
irradiation effects for any reasonable system ages. Since we are not performing evolutionary 
calculations (Burrows et al. 2000; Burrows, Sudarsky, & Hubbard 2003 (BSH); Baraffe et 
al. 2003) in this paper, the system ages are not germane to the problem at hand, which is 
reproducing the observed planet-star flux-density ratios in the Spitzer bands. These ratios 
depend only on the stellar spectrum, the orbital distance, and the planet’s mass and radius, 
but not directly on age. 

The KOV stellar primary of the transiting extrasolar giant planet TrES-1 is 157 ± 6 
parsecs distant, has a Tefjof 5214 ± 23 K, a metallicity near solar {[Fe/H] ~ 0), a radius 
(R*) of 0.83 ±0.03 Rq, a mass (M*) of 0.87±0.05 M©, and a bolometric luminosity (L*) near 
half solar (Laughlin et al. 2005). The planet’s orbital and physical parameters are a semi¬ 
major axis (a) of 0.0393 AU, a period (P) of 3.030 days, a planet mass (Mp) of 0.729 ±0.036 
Mj ^ (Laughlin et al. 2005) or 0.76 ± 0.05 Mj (Sozzetti et al. 2004), and a transit radius 
(BSH; Burrows et al. 2004) of approximately 1.08 ± 0.05 Rj ^ (Alonso et al. 2004; Sozzetti 
et al. 2004; Laughlin et al. 2005). 

The corresponding quantities for the F8V/G0V star HD209458 are 47.3 parsecs (Perry¬ 
man 1997), Teff -6000 K, [Fe/H] ~ 0, R* = 1.2±0.1 R©, M,= 1.1±0.1 M©, and L, = 1.6 L© 
(Henry et al. 2000; Charbonneau et al. 2000; Brown et al. 2001). The planet’s parameters are 
a = 0.0468 AU, P = 3.524738 days, Mp = 0.69 ±0.02 Mj , and Rp = 1.32 to 1.40 Rj (Henry 
et al. 2000; Charbonneau et al. 2000; Mazeh et al. 2000; Brown et al. 2001; Cody and Sasselov 
2002; Fortney et al. 2003; Laughlin et al. 2005), with the lower value slightly more favored. 
Both TrES-1 and HD209458b have orbital inclinations near 90° and eccentricities near 0. 

We see that the flux at the surface of TrES-1 is less than half that at the surface of 
HD209458b. Eurthermore, the larger radius of HD209458b implies that it intercepts more 
than two times the radiative power. The former implies that the atmospheric temperature 
of HD209458b is likely to be higher than that of TrES-1, while the latter implies that 
HD209458b is bolometrically brighter. Nevertheless, as we show, their theoretical planet- 
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star ratios are not expected to be very different. 

Prior to the recent measurements of the secondary eclipses that motivate this paper, 
Charbonneau et ah (2002) had uncovered evidence for the presence of sodium in the at¬ 
mosphere of HD209458b from the different transit radii in and out of the Na-D line (see 
also Fortney et al. 2003). Similarly, Vidal-Madjar et ah (2003) had seen evidence for 
atomic hydrogen in a stellar-flux-induced planetary wind from HD209458b. The transit dip 
in Lyman-a is of such a magnitude (~15%) that the planetary material clearly extends be¬ 
yond the Roche lobe, indicating planetary mass loss (Burrows & Lunine 1995; Lecavalier des 
Etangs et al. 2004). Other than these data, there had been no determinations to date of the 
composition of an extrasolar planet. 

The new secondary eclipse data for TrES-1 (Charbonneau et al. 2005) and HD209458b 
(Deming et al. 2005) hint that this situation may be changing. TrES-1 shows eclipse 
depths (planet-star flux-density ratios) in the Spitzer/IRAC band centered at 4.5 /xmof 
0.00066±0.00013 and in the Spitzer/IRAC band centered at 8.0 /rmof 0.00225±0.00036 
(Charbonneau et al. 2005). HD209458b shows a corresponding ratio of 0.0026±0.00045 in 
the Spitzer/MIPS band centered near 24 pm (Deming et al. 2005). These numbers are actu¬ 
ally the ratios of detected electrons, an approximate substitute for the ratio of average flux 
densities in a given band. For flux density comparisons, one relies on flux calibrations that 
may not yet be robust, particularly given the signihcant differences between the spectra of a 
close-in EGP and a calibration star (e.g., Vega). Given this, in §4 we also provide theoretical 
bandpass-averaged detected-electron ratios. 

Since the data from Charbonneau et al. (2005) and Deming et al. (2005) are but three 
of the ten potential data points (hve bands x two nearby transiting planets) we can expect 
using Spitzer, very-low-resolution spectra (but “spectra” nevertheless) of EGPs are antici¬ 
pated soon that will provide compositional and atmospheric information of an unprecedented 
character. 


3. Numerical Techniques, Databases, and Assumptions 

The numerical tools we employ to derive the close-in planet’s spectrum during secondary 
eclipse are described in BSH and Burrows, Sudarsky, & Hubeny (2004), to which we refer 
the reader for further details. The spectral/atmosphere code COOLTLUSTY (Hubeny 1988; 
Hubeny and Lanz 1995; Sudarsky, Burrows, & Hubeny 2003) handles the effects of external 
irradiation using a hrst-order variant of the DFE (Discontinuous-Finite-Element) method. 
The incident flux is isotropically spread over the hemisphere of a given planar patch of 
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planetary “surface.” As described in BSH, to account in approximate fashion for the variation 
in incident flux with latitude when using a planar atmosphere code, as well as the possible 
day-night differences, we introduce the flux parameter /. A value of / = 0.5 assumes 
that there is little sharing of heat between the day and night sides of the EGP. A value 
of / = 0.25 assumes in the calculation of the T/P prohle that the heat from irradiation is 
uniformly distributed by efficient winds over the entire sphere and that the infrared emissions 
are isotropic. For this study of the planet-star flux-density ratios of close-in EGPs, we use a 
hducial value for / of 0.25, but we return to the issue of the proper choice of / in §5. 

The emergent planetary spectra are calculated at 5000 wavelength points distributed 
logarithmically from 0.3 to 300 microns. The temperature-pressure prohle in hux equilibirum 
is derived from ~10“® to ~10^ bars. The spectral models for the stellar primaries (given a 
specihc spectral subtype) come from Kurucz (1994). The molecular and atomic opacities are 
taken from the opacity library described in Burrows et ah (2001) and equilibrium composi¬ 
tions are derived using the updated thermochemical database of Burrows & Sharp (1999). 
Though the atmosphere is close to an ideal gas, the H/He equation of state of Saumon, 
Ghabrier, & Van Horn (1995) that incorporates non-ideal effects is used. 


4. Comparison of Theoretical Planet-Star Flux Ratios 

Figure 2 depicts our theoretical planet-star hux-density ratios versus wavelength in the 
near- and mid-infrared for TrES-1 (magenta) and HD209458b (green). To derive these curves 
we have used the physical data for the planets and their primaries described in §2. Our 
baseline planet models have solar metallicity. The phase-averaged (Sudarsky, Burrows, & 
Hubeny 2003) flux ratios, but for / = 0.25 (§3; Sudarsky, Burrows, & Hubeny 2003; Burrows, 
Sudarsky, & Hubeny 2004), are given and have not been shifted in any way. Superposed are 
the new data at 24 /xmfor HD209458b (green) and at 4.5 /xmand 8.0 yumfor TrES-1 (gold). 
The vertical error bars are the quoted 1-a ranges and the horizontal bars indicate the widths 
of the corresponding IRAG and MIPS bands. Also included in yellow are the four theoretical 
IRAG band-averaged fluxes for the TrES-1 model and in blue for the HD209458b model. 
These band averages are derived using the published transmission functions divided by 
frequency to obtain the theoretical ratio of detected electrons. Since the close-in EGP (this 
work) and stellar (Kurucz 1994) spectra are so flat at 24 /rm, a transmission band-averaged 
point is not shown (or needed) there. 

Varying Mp, Rp , and R* within the error bars alters the resulting planet-star flux-density 
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ratios only slightly. Similarly, and perhaps surprisingly, adding Fe and forsterite clouds does 
not shift the predictions in the Spitzer bands by an appreciable amount. Moreover, despite 
the more than a factor of two difference in the stellar flux at the planet, the predictions 
for the planet-star ratios for two such disparate close-in EGPs as TrES-1 and HD209458b 
are not very different. For TrES-1, raising or lowering the metallicity by a factor of three 
changes the flux ratios by less than ~10%. Interestingly, however, changes in / and, by 
inference, day-night atmospheric differences and phase function effects can result in 25% to 
50% deviations in the flux ratios of which one should take note. 

Comparing our baseline models with data on Fig. 2, we can deduce several interesting 
things. First, the 24-/imdata point is close to the predicted value, though for all three data 
points the theory slightly underestimates the data by a factor of ~1.5-1.8, with the largest 
discrepancy being for the TrES-1 band at 8.0 pm. For HD209458b at 24 pm this deviation 
is only ~ la. The high absolute values of the flux ratios imply that the close-in EGP 
atmospheres are indeed at high temperature, predicted for TrES-1 and HD209459b to be 
~1500 K and ~1600 K, respectively, at a Rosseland depth of ~1. Both atmospheres clearly 
span the temperature range 1000-2000 K. For TrES-1 at 8 pm, theory yields a brightness 
temperature (temperature at t\ = 2/3) near 800-900 K, slightly lower than the ~1100 K 
crudely inferred from the data. At 4.5 pm, the theoretical brightness temperature of TrES-1 
is ~750-900 K (using / = 0.25), again slightly lower than the data might imply. However, 
care must be taken in estimating temperatures of any sort, and we will not, due to the 
hazards of extracting an “effective” or “equilibrium” temperature from these data, say much 
more about them. We do note, however, that given the radius of HD209458b and its general 
flux level, we expect that its bolometric luminosity is above 2 x 10“^ Lq. This is approaching 
the luminosity of a star ~100x the mass of HD209458b at the edge of the hydrogen-burning 
main sequence (Burrows et al. 2001). 

In Fig. 2, there is a hint of the presence of H 2 O, since it is expected to suppress flux 
between 4 pm and 10 pm . This is shortward of the predicted 10-pm peak in planet-star flux- 
density ratio, which is due to water’s relative abundance and the strength of its absorption 
bands in that wavelength range. Figure 3 demonstrates this by comparing H 2 O, GH4, and 
GO opacities per molecule from 0.5 to 10 microns. Without H 2 O, the fluxes in the IRAG 
bands would be much higher than the fluxes in the mid-infrared. Hence, a comparision of 
the TrES-1 and HD209458b data at 4.5/8.0 pm and 24 pm suggests, but does not prove, the 
presence of water. As Fig. 3 implies, seeing the expected slope between the 5.8-pmand 
8.0-pmbands and the rise from 4.5 pm to 3.6 pm would be more revealing in this regard and 
is a prediction of our theory. Furthermore, the relative strength of 24-pmMIPS flux ratio in 
comparison with the 3.6-pm, 4.5-pm, and 5.8-pm IRAG channel ratios is another prediction 
of the models, as is the closeness of the 8.0-pmand 24-pm ratios. The latter seems borne 
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out by the data, though these data are for two different objects. 

Nevertheless, the models have difficulty htting the depth of the 4.5-/rm feature in TrES-1. 
This feature coincides with the strong CO absorption predicted to be a signature of hot EGP 
atmospheres (Fig. 3; Sudarsky, Burrows, & Hubeny 2003; Burrows , Sudarsky, & Hubeny 
2004; Burrows 2005), but is shallower than expected and ~2-(T discrepant. The data are 
in fact the band-averaged flux-density ratios of the detected electrons. As such, the larger 
fluxes on either side of the trough theoretically centered close to 4.5 /rm contribute planet 
flux to the detected band. As a result, the yellow dot that represents this integrated band 
contribution is a weak function of CO abundance. In fact, a CO abundance 100 x larger than 
expected in chemical equilibrium lowers this flux ratio at 4.5 /xmby only ~25%. Therefore, 
while the 4.5-/imdata point for TrES-1 implies that CO has been detected, the exact £t 
is problematic. More data and further attention to calibration are called for. However, as 
we have indicated and discuss in §5, the contrast between theory and measurement at all 
data points may be a signature of day-night infrared flux asymmetry. The close-in planets 
may not be radiating heat energy isotropically (as is assumed when using / = 0.25), a 
not-unexpected result (Guillot & Showman 2002). 


5. Discussion 

We predict that both the 3.6-/imand the 5.8-/imflux ratios will be higher than the 
4.5-/imratio by at least 50% and that the pattern of the yellow and blue dots on Fig. 2 
will be realized. We also predict that the 24-/imflux ratio for TrES-1 will be similar to 
that seen for HD209458b (Doming et al. 2005). However, while the close correspondence of 
the measured and theoretical fluxes depicted on Fig. 2 is striking, it is not perfect. What 
could explain the differences? One major uncertainty is the day-night atmospheric profile 
difference. HD209458b and TrES-1 are close enough to their primaries to be in synchronous 
rotation. Therefore, they show the same hemisphere to the star at all times. It is the 
zonal winds, atmospheric circulation currents, and jet streams (Menou et al.2002; Guillot & 
Showman 2002; Showman & Guillot 2002; Cho et al. 2003; Burkert et al. 2005) that advect 
heat from the day to the night sides, thereby affecting the atmospheric temperature structure 
as a function of longitude. How much of the stellar radiation goes to heating the day side 
(visible just before and after the secondary eclipse) and how much is transported by mass 
motion away from the day side to heat the night side? This question remains unresolved, but 
directly impinges upon the planet-star flux-density ratios in the infrared measured during 
secondary eclipses. 

The / factor we use for our fiducial model (0.25) is tailored to distribute heat over 



the entire planet, on both the day and the night sides. The three data points depicted 
in Fig. 2 are all factors of ~1.5-1.8 above the / = 0.25 theory curves. This implies that 
the planet’s radiation is predominantly radiated by the day side and that its emissions 
are forward-beamed. This is not entirely unexpected in the optical, but seems to be the 
case in the mid-infrared as well. The implication of the new secondary eclipse data at 24 
/imand 8 /rmmay be that we are seeing indirect signs of an asymmetry in the day-night 
heating and temperature profiles, with the day side hotter than the night side by at least 
500 K. This estimate is based on the mid-infrared “excesses” seen in Fig. 2, on the possible 
model variations, and on the measurement errors. However, the conhrmation of such a 
conclusion awaits more detailed multi-dimensional general circulation models with correct 
transport, next-generation models of the wavelength-dependent phase functions of close-in 
EGPs (Sudarsky, Burrows, Hubeny, & Li 2005), and, most importantly, additional data. 

It has long been suggested, and recently articulated (Cooper & Showman 2005), that 
due to winds the sub-stellar point of an irradiated close-in EGP may not be the hottest spot. 
Advection would introduce a lag even for circular orbits between the planet’s ephemeris and 
its light curve (the light curve would lead). Cooper & Showman (2005) estimate that the lead 
could be as much as 60° and could amount to a 20% brightness shift in the value at superior 
conjunction. A 20% decrement is not enough to close the modest apparent gap between 
our theory and the 24-/rmdata for HD209458b, and is not of the correct sign. However, the 
concept of a shift in the light curve deserves further scrutiny. 

Charbonneau et al. (2005) estimate a Bond albedo (A) for TrES-1 of 0.31±0.14. How¬ 
ever, one should be very cautious using these new data to infer temperatures and reflection 
coefficients. Not only is the reradiation not expected to be isotropic off the planet, but the 
atmospheres are not black bodies. While one can distinguish hot atmospheres (1000-2000 
K) from cooler atmospheres (500-1000 K), the data and theory are not yet adequate to allow 
these TrES-1 data to strongly constrain A. However, if A were in the 0.31±0.14 range, this 
would imply that there is a cloud of non-trivial optical depth in the upper layers of TrES- 
1, putting it into the “Class V” category of EGPs, rather than the “Class IV” category 
(Sudarsky, Burrows, & Pinto 2000). The latter, due to strong absorption bands and little 
Rayleigh scattering, have very low Bond albedos (below 0.05). Our models for TrES-1 would 
favor the low-albedo Class IV category for TrES-1 and the higher-albedo Class V category 
for HD209458b, but we feel it is premature to conclude anything dehnitive about albedos at 
this stage (other than that they can not be very high and that the close-in EGPs can not be 
highly reflective). 

In this paper, we have calculated planet-star flux-density ratios versus wavelength, fo¬ 
cussing on the near- and mid-infrared out to 25 microns and the irradiated close-in extrasolar 
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giant planets TrES-1 and HD209458b, and have compared our theory with the recent sec¬ 
ondary eclipse data from Charbonneau et ah (2005) and Deming et al. (2005). We have 
inferred the presence of carbon monoxide, and perhaps water, in the atmosphere of TrES-1 
and have determined that both atmospheres are hot (^1000-2000 K). We have explored the 
effects of varying the metallicity, Mp, Rp, and R*, the latter three within the stated error 
bars, and find our predictions for the planet-star flux-density ratios to be robust. However, 
more work is required to understand the apparent shallowness of the TrES-1 4.5-/im trough. 
(Non-LTE and non-equilibrium effects in the chemistry may be worth exploring, as well as 
the consequences of possible stratospheric temperature inversions (Smith & Hunten 1990).) 
However, we suggest that the data (particularly at 8 and 24 microns) indicate we are begin¬ 
ning to constrain the degree of anisotropy in the temperature prohle of a close-in EGP (the 
day-night contrast) and in the angular dependence of its emission. Our preliminary conclu¬ 
sion is that the slight systematic differences seen in Fig. 2 between the phase-averaged theory 
and all three new measurements may be explained in part by an infrared-brighter, hotter 
day side. This raises the intriguing possibility that additional and more precise secondary 
eclipse data for these transiting EGPs could shed light on their global meteorology. Be that 
as it may, these secondary eclipse data are opening a new chapter in the accelerating study 
of extrasolar planets and emphasize that knowledge of a unique character is generated when 
new capabilities emerge. 
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Mass (Mj) 


Fig. 1.— Transit radii (in units of Jupiter’s radius) with error bars versus planet mass with 
error bars (in units of Jupiter’s mass) for the seven EGPs currently seen to transit their 
primaries. The positions of Jupiter and Saturn themselves are included for comparison. 
Note that HD209458b is the largest transiting EGP (Mazeh et ah 2000; Brown et ah 2001; 
Gody and Sasselov 2002) and, as such, is an outlier, though most of these irradiated EGPs 
are clearly larger than Jupiter. An extended transit radius is a known consequence of stellar 
irradiation (Guillot et ah 1996; Burrows et ah 2000; Burrows, Sudarsky, & Hubbard 2003; 
Baraffe et ah 2003; Ghabrier et al. 2004; Burrows et ah 2004). 
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X (microns) 


Fig. 2.— The logarithm base ten of the planet-to-star flux-density ratio as a function of 
wavelength (A, in microns) for our baseline models of TrES-1 and HD209458b (for / = 
0.25). The model for TrES-1 is purple and that for HD209458b is green. Superposed are 
the secondary eclipse data: the gold dots with corresponding error bars are the TrES-1 
Spitzer/IRAC data from Charbonneau et ah (2005), while the green dot with error bars is 
the HD209458b Spitzer/MIPS 24-/rmdatum from Deming et al. (2005). Also included are 
the band-averaged detected-electron/“flux” ratios for the TrES-1 (yellow) and HD209458b 
(blue) models in the four IRAC bands. Note that coincidently the blue dot at 4.5 pm overlaps 
the gold TRES-1 data point. The position of the strong CO absorption feature at ~4.67 
pm is indicated and clearly coincides with the ~4.5 pm IRAC band flux. See text for a 
discussion and details. 
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Fig. 3.— The logarithm base ten of the absorption cross section per molecule (in cm^) versus 
wavelength (in microns) from 0.5 /imto 10 /rmfor water (H 2 O), carbon monoxide (CO), and 
methane (CH4). These numbers are not weighted by abundance. A comparison with Fig. 2 
clearly shows that CO is probably present in TrES-1 and that water helps shape the planet- 
star flux-density ratio in the 4 to 9 micron spectral region. Shown also are the positions of 
various standard photometric bands {Z, J, H, K, M). This hgure, in conjunction with Fig. 2, 
demonstrates that if CH4 is present in abundance in either HD209458b or TrES-1 then the 
3.6-/imIRAC band will test this. However, our expectation for these close-in EGPs is that 
CH4 will not be in evidence. When data for all four IRAC bands are available, we should 
be able to verify dehnitively the presence of water vapor in the atmosphere of an extrasolar 
giant planet, though by the absolute level at 8.0 /imof the TrES-1 flux in Fig. 2 the presence 
of water already seems likely. See text for further discussion. 






















